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Abstract 


We  present  here  O(logn)  time  parallel  algorithms,  on  the 
CREW  PRAM  model  of  parallel  computation,  for  determining  all 
the  arcs  and  their  sequence  in  the  boundary  of  the  intersection 
or  the  union  of  n  circular  discs  in  the  plane.  Here,  the  inter- 
section of  n  circular  discs  is  the  part  of  the  plane  which  be- 

3 
longs  to  all  of  them.  Our  techniques  use  0(n  )  processors.  Our 

algorithms  take  O(loglog  n)  time  in  the  CRCW  PRAM  model  of  pa- 
rallel computation. 

We  also  present  matching  lower  bounds  for  the  parallel  ti- 
me of  the  intersection  and  union  problems.  Our  algorithms  are 
based  on  an  interesting  characterization  of  the  boundary  of 
the  intersection  and  the  union.  We  also  show  that  the  boundary 
of  the  intersection  of  n  circular  discs  in  the  plane,  consists 
of  at  most  2n-2  arcs,  for  any  n'^2 . 
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Introduction 


The  problem  of  computing  the  intersection  or  union  of  a 
set  of  planar  objects  has  been  given  some  attention  in  the  past 
(see  e.g.  [Bentley,  Ottman,  79],  [Sharir,  83]).  When  the  objects 
are  circular  discs  or  boxes  then  the  problem  relates  to  that  of 
quick  computation  of  intersection  or  union  of  relations  in  da- 
tabase applications. 

Related  problems,  such  as  estimating  the  area  of  the  union 

of  many  discs  in  the  plane,  were  first  posed  by  [Shamos,  78]. 

Monte-Carlo  techniques  for  computing  the  volume  of  the  union 

of  n  spheres  in  k  dimensions  were  given  in  [Spirakis,  85].  They 

run  in  time  0(nk)  and  are  based  on  a  method  developed   by  [Karp, 

Luby,  83]  for  estimation  of  the  failure  probability  of  an  n- 

component  system.  Also,  [Sharir,  83]  shows  how  to  construct  the 

2 
Voronoi  diagram  of  a  set  of  n  circular  bodies,  in  time  O(nlog  n) 

and  this  can  be  used  for  detecting  intersections  between  the 
circular  bodies.  [Brown,  79],  shows  O(nlogn)  sequential  algo- 
rithms, based  on  inversion  transforms,  which  compute  the  (boun- 
dam  of)  the  intersection  and  union  of  n  circles  in  the  plane 
(see  also  [Mehlhorn  ,  198'^]). 

We  consider  here  the  problem  of  computing  all  the  arcs  and 
their  sequence  in  the  boundary  of  the  intersection  (or  the  union 
of  n  circular  discs  in  the  plane.  We  first  provide  a  nice  cha- 
racterization of  those  boundaries.  This  leads  to  an  O(logn) 
parallel  time  algorithm,  running  on  the  CREW  PRAM  model  of  pa- 
rallel computation.  Our  algorithm  takes  O(loglogn)  time  on  a 

3 
CRCW  PRAM.  It  needs  n   processors.  Since  O(nlogn)  sequential 

algorithms  for  computing  the  intersection  or  union  are  known, 
our  bound  on  the  number  of  processors  is  by  no  means  the  best 
possible,  (see  [Brown,  79]).  We  also  show  that  the  boundary  of 
the  intersection  of  n  circular  discs  in  the  plane,  consists  of 
at  most  2n-2  arcs,  for  any  n^2 .  To  the  best  of  our  knowledge, 
this  is  the  first  set  of  results  concerning  the  parallel  compu- 
tation of  the  boundaries  of  union  and  intersection  of  many  discs 
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Notation  and  definitions 


se 


A  circular  disk  K,  of  radious  r  and  center  (x,y)  is  the 

t  ^( X ' , y ' ) : ( x-x ' )  +(y-y')  4 r  y    Its  interior,  K,  is  the  set 

t2        2   2  1 
(x',y'):(x-x')  +(y-y')<r  jand  its  circumference,  C,  is  the 

set  ^  (x'  ,y • ) :  (x-x'  )^+(y-y ' )^  =  r^J .  It  follows  that  K  =  KUC  and  that 

K«C  =  0. 

All  arcs  in  C  are  directed  counterclockwise  (in  this  paper), 
A  closed  arc  (which  includes  the  end-points  a,b)  is  denoted 
[a,b].  Parentheses  are  used  for  open  and  semi-open  arcs.  So,  the 
complement  of  [a,b],  is  C-[ a , b ] = ( b , a ) =[ a , b ]  .  An  one  point  arc 
[a,a]=  a   has  (a, a)  as  its  complement,  hence  C=  ( a ^^ U{a,a)=[a,a) 
=  (a ,a]  . 

The  ordering  of  points  on  the  circumference  C  of  K,  in  the 

counterclockwise  direction,  depends  on  the  point  of  C  used  as 

an  origin.  In  Figure  1,  q  preceeds  r  if  p  is  the  origin  while 
q  follows  r  if  p'  is  the  origin. 

Definition   We  say  that  point  x  in  C  preceeds  point  y  in  C  with 
respect  to  p  in  C  iff  x  6[p,y].  We  denote  this  by  p:^x,y> 

More  generally,  let  L  be  a  finite  list  of  points  in  C  which 
may  include  p  and  may  also  include  multiple  points. 

Definition   We  say  that  L  is  ordered  with  respect  to  p  G  C,  and 


we  denote  this  by  L=p<x  , 


,Xj^>  iff  p:<x^,x  >  for  i^j  and 


i  ,  j^k. 

Definition   Let  L  be  a  finite  list  of  points  in  C.  Consider  L 
as  a  set  i.e.  drop  repeated  points.  We  denote  by  Right  (x;Ll 
the  point  r  of  C  with  the  property 

[x,r]f\L  =  [r\ 
Also,  we  denote  by  Left  (x;L)  the  point  h  of  C  such  that 

[h,x]OL={h\ 

Note   If  xGL  then  Left  (x;L)=Right  (x;L)=x.  If  L= x <x , , x, , . . . , x  S 
and  all  x^'s  are  distinct,  then  Right(x;L)=x   and  Lef  t  (  x  ;  L  )  =Xj^ 


Now,  we  are  given  n  discs,  K.,  i=1,...,n  (with  circum- 
ferences C.,  i=1,...,n).  Without  loss  of  generality  assume 

K.^K.  for  i^fj. 
1   J        ^ 

Definition 

n 

The  intersection  of  K,,...,K   is  the  set  1=0  K. 

1     '  n  .^^1 

n 

The  union  of  K",  ,  .  .  .  ,K   is  the  set  1=  "IT  K ,■  . 

In  •  ,  •'• 

1  =  1 

Definition 

Let  us  fix  i.  If  K.f^K.  i    0  for  some  j,  then  consider  the 
_         ^   J  _ 

unique  arc  C.r>K.;^0,  which  is  the  part  of  C.  that  lies  in  K.. 

Let  b.  .  be  the  beginning  and  e.  .  the  end  of  this  arc  (ordered 

counterclockwise  on  C.). 

1   _ 

Thus   [b  .  .  ,e  .  .  ]  =  C.f\K  .  £   c  . 

Clearly  (Fig.  2)  b..=e.  ..  The  first  index  is  always  meant 
with  respect  to  the  circumference  which  the  arc  belongs  to. 

Thus,  [b..,e..]  is  an  arc  in  C.,  while  re..,b..1  is  an 
arc  in  C. ,  both  having  the  same  endpoints. 

The  intersection 


A.=  ry  (C.AK  .)=  n  [b.  .,e   ]  c  C 

1  j^i   1   J   j^i   ij   ij   -   1 

is  the  part  of  C.  which  lies  within  all  K.,  -[ii. 

1  J.J 

empty  or  consist  of  one  or  more  arcs. 


(») 


A.  could  be 


Similarly,  the  union 


V.  =  \r  (C.nK  .) 
'  j^i   ^   J 


is  such  that 


its  complement,  with  respect  to  C.,  i.e. 

V%  n  (C  AK  )^=  n  (e.  .,b.  .) 

^  j'^i   '   J    j^i   iJ   iJ 

is  the  part  of  C.  that  lies  outside  all  K.,  j/i 


*  « 


i=1 


A.  Characterization  of  the  boundaries 

c 

Intuition  says  that  A-  (respectively  V.)  must  be  the  contri- 

1       n  _        ^  n  _ 

bution  of  C.  to  the  boundary  of  n  K.  {respectively  of  ^S   K.). 

^  i  =  1  ^ 

The  following  theorem  proves  this: 

Theorem  A . 1  n 

(a)  The  boundary  of  the  intersection  I  is  "U*  A^  • 

n   ^   i  =  1 

(b)  The  boundary  of  the  union  I  is  "VT  V .  . 


i=1 


Proof 


We  show  only  part  (a)  here.  Part  (b)  is  similar. 

n        n 
It  is  enough  to  show  that  I=(  f\    K.  )U(  V  A.  ) 


i  =  1 
n        n 
i)  Let  I'  be  the  set  {   (\    K.JUlVA.) 


i=l 


i=l  ^ 


i  =  l 


Let  xei.  This  means  x  €  (  K  .U  C  .  )  ^^  i  =  1  ,  .  .  .  ,  n  . 

11 

If  xSK .  for  all  i=1,...,n  then  clearly  x  6  I' 
1 

Now  consider  x  €  C   for  some  m. 

m 

Since  also  xGI,  we  get  that 


x  e  (I  f^c  ) 

m 


I.e. 


6(C  n  (nK,))=^    x6A   r^x6l' 
m    i-  T  1  m  ^ 


(ii)  Conversely,  assume  xGI  • 

n        n 
I.e.    xe  (nK.)U(VA 


i=1 


i=l 


So,  either  x  £(  H  k  )  or  x  €  A   for  some  m 
i=l  i  " 


n 

If  X  6(  n   K 

i  =  l  ■ 


n       n 

Then  x  6  I,  since  0  K  S  H  K 

i=1  '   i=1  i 


If  x  eA   for  some  m  then  x6  il 


m 


C  OK  .  ) 


But       C  ft  K  .  £.  K  .         and    C     £  K 
ra      J         J  mm 


So       X    6     0    K  .  =  1 
J=1     J 


QED 


Theorem  A.I  is  useful  to  a  characterization  of  the  sets 
c 


A.  ,  V 

1    1 

Definition 


We  define  the  lists  B.,E.,L.  of  points  of  C.  to  be 

1  '  1   1     ^  1 

B.  :  =  (b.  J,  j/i,  j=1  ,  .  .  .  ,n) 


E^:  =  (e.  ^,  j/i,  j=1  , 


L. :=(b.  . 
1     ij 


IJ 


.,n) 
,  j;^i  ,  j  =  1  ,  .  .  .  ,n 


All  lists  may  contain  multiple  points. 

Note   If  C.ftK.  =  0  for  all  i  ;^  i  then  all  lists  are  empty. 
1   J  ^ 

If  C.rNK"  =0  for  some  m,  then  A.=0,  v'^  =  C.. 

1   m  111 

Definition 

P(i)  is  defined  to  be  the  set  of  indices  j  satisfying 


b.  .:<e.,  ,b.,>  for  all  k^^i.i 
1 J    ik   ik'  '  -^ 

Similarly, 


:with  b . ,  ^b  .  . 
ik   1 J 


e{i)  is  defined  to  be  the  set  of  indices  j  satisfying 

e..:<b.,,e.,S  for  all  k;ii,i   (with  e.,;i^e..) 
ij    ikik'^  '  -^  ikij 

The  following  two  theorems  compose  the  main  result  of  our 

paper  and  lead  to  a  fast  parallel  algorithm  for  computing  the 

boundaries  of  I  and  I: 

Theorem  A  .  2 
(a)  A. = 


i=  V   [b.  .  .Right  (b.  .  ;E.  )  ] 
i=  "^    [e.  .  .Right  (e.  .;B.  )  ] 


Proof 


Part  a 


T •  A^S  V"    [b.  ..Right  (b.  .,E.  )] 
^   je(3(i)   iJ    ^      ij'  1 


Assume  x  6A..  Then  there  exists  a  j  such  that  b^  =Lef t ( x ; B^ 
Let  e^p  be  the  Right  (b.  .;E^). 
We  will  show  that  x  e[b.  -.e.  ] 
We  consider  two  cases 

Case  1   If  b..=x  then  clearly  x  €[b..,e   ] 
ij  ij   ir 

Case  2   b  .ix.    Then  x;i^b.,  for  all  k,  by  the  definition  of  Left. 

Also,  (for  same  reason)  b.   6  (x,b.  .J. 

ir      '  ij 

Now,  since  x  6  A.  we  get 


X  e  [b.  , e .  ] 

ir   ir 


So,  we  have 
and 


b  .  :  <  X  ,  e  .  > 

ir       ir 

b.  :  <  b.  .  ,x> 
ir  ^   1 J ' 


EQ1  ) 

by  EQ1  ) 
by  Left 


Hence  ,    b.  :<'b..,x,e.  "> 
ir^ijir' 

(by  transitivity  of  "less  than") 

I.e.     xe[b..,e.  ]. 

1 J   ir 


We  still  have  to  show  that  j  6  P(i). 

Assume  the  opposite,  i.e.  that 

3l<;^i,i  with  b  .  ,  ih  .  .    such  that  b.,6(b..,e.,) 
'-*  ~^  ik   1 J  ik    1 J   ik 

Remembering  how  b.  .  was  defined,  we  get   x:^b.,  ,b.  .S. 

But  then,  x  %    [b.,  ,e.,  ]   =>   x  8  A.   which  is  a  condradic tion . 

ik   ik  1 


So,  we  showed  t 


hat  AC  V*    [b.  .,Right(b.  .;E.  )  ] 


□ 


2-   'U"    [b.  .,Right(b.  .;E.  )  ]£  A. 
j6(3(i)   iJ    ^     ij'  1      1 

Since  j  6  p(i)  we  get  (by  def.  of  p) 

[b.  .  ,e  .  ]S  [b.  ,e.  ]   V   m;ii         (EQ3 
1 J   im      im   im 

Also,  by  definition  of  Right: 

[b.  .,Right(b.  .;e.  )  ]£  [b.  .,e.  ]   V^  m   (EQ^) 
1 J         1 J  '  1  ■'   ^  ij  '  im' 

by  (EQ3)  and  (EQA)  we  get 

[b.  .  , Right (b.  .  ;E.  ) ] S  [b.  ,e.  ]   V  m 
1 J    ^     1 J   1       im   im 
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I.e.   [b.  .,Right(b.  .',  E.)1Ca. 
and  this  holds  for  all  j6P(i) 
Hence    -y      [ b.  , Right ( b^  ; E^ ) ] 
j6p(i) 


£  A 


□ 


(1)  and  (2)  prove  part  a  of  the  Theorem. 
Part  b   Proof  in  full  paper 

Theorem  A. 2  alone  would  not  be  so  useful  if  the  intervals 
were  overlapping.  Fortunately,  the  next  theorem  shows  that  this 
is  not  the  case. 

Theorem  A . 3 

(a)  Let  A .  =[b.  , Right (b .  ;E .  )  ] 

im    im         im   1 

for  m  =  j,k  6  (3(i) 

Then  either  A..  =  A.,   or  A  .  .CX  h  .     =<Z 
ij   Ik        ij    Ik 

(b)  Let  r .  =[e .  , Right (e  .  ;B.  )  ] 

im    im    °     im   1 

for  m=j,k  e  e{i) 

Then  either  T.  .=r.,   or 
ij   Ik 

For  proof  see  Appendix  A1 


r.  .n  r.,  =0 

IJ    Ik 


5.  O(logn)  parallel  time  Algorithms  for  Intersection  and  Union. 

3 
We  use  a  total  of  0(n  )  processors,  allocated  in  groups  of 

2  — 
n   per  cycle.  The  processors  of  cycle  K.  are  allocated  as  fol- 
lows: A  distinct  group,  G.  .,  of  n  processors,  one  for  each  pair 
K.  ,  K.     {j^D. 

1.  Initialization 

Initially,  we  let  a  specific  processor  of  each  group  G.  . 
to  compute  (in  parallel  for  each  i  and  i)  the  two  points  b.  ., 

e    of  the  intersection  of  the  two  circumferences  C.,C..  The 

1 J  1   J 

b   ,e    are  computed  by  finding  their  angles  (with  respect  to 
a  direction  9=0   independent  of  the  discs)as  seen  from  the  cen- 
ter of  disk   K..  Let  these  angles  be  9(b. .),  S(e.  .)  as  in  Figu- 
1  ij  '     ij 

re  3  . 


This  computation  takes  a  constant  number  of  parallel  steps. 
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2  .  Computing  P(  i  )  ,   £(  i  ) 

Initially,  a  distinct  processor  of  the  group  G^  initiali- 
zes a  shared  variable,  check  , to  true.  Then,  each  group  G^  of 
processors  checks  in  parallel  whether  j  belongs  to  p(i). 

That  is,  each  processor  p  (k=l,...,n)  of  the  group  G^., 
checks  in  parallel  whether  ^^  .  :<e^|^ ,  b^l^>  for  b^  /b^^^.  Processors 
p,^  record  violations  of  b  .  .:<e  .  ^^  ,h  .  ^    by  setting  a  local  bit  B^ 


ij 


ik'  ik' 


to  false.  (If  p,  records  no  violation,  then  B,  is  set  to  true), 
k  k 

Finally,  the  logical  AND  of  the  B   is  computed  in  check   . 

K  1  J 

This  takes  O(logn)  time,  in  the  CREW  PRAM  model,  and  0(1)  time  in 
the  CRCW-PRAM.   A  similar  procedure  is  followed  for  the  case  of 
computing  e(i).  Processors  p   now  record  violations  of  ^i-i=^bj_^;' 
e  "^  and  again  a  logical  AND  has  to  be  performed  in  parallel. 

3.  Finding  A.  (and  V*^) 
°   1         1 


Each  arroup  G.  .  reads  check.  ..  Conditioned  on  check.  .  =  true, 

G.  .  constructs  the  arc  of  C.,  with  start  b.  .,  which  contributes 

to  A-  (with  start  e.  .,  if  it  contributes  to  V.)  as  follows: 
1  1 J  1 

Each  processor  p,  in  G.  .  computes  the  angles 
^  k      1 J 

a.,  =3(b.,  )-3(b.  .) 
Ik     ik      1 J 

T.,  =9(e.,  )-9(b.  .) 
ik     Ik      1 J 

Then  the  group  G.  .  comoutes  the  x.  .=min)T-,  ,    k;^i\   and 
^  .  ^   ij     •  ij     <•  Ik'     ■» 

y.  .=minfa.,  ,  k/iL  This  takes  O(logn)  time  in  the  CREW  PRAM  mo- 
•'ij     Ilk      ■»  ° 

del.  The  contribution  to  A.  is  [b.  .,x.  .]  and  to  V.  it  is  [e.  .,y.  . 

1       ij '  ij  1  1 J   1 J 

(In  the  CRCW  PRAM  model,  step  3  would  take  O(loglogn)  time). 

Note 

The  algorithm  just  presented  constructs  the  contour  surroun- 
ding the  intersection  I  (if  li<Z)    and  the  union  E  of  the  n  circu- 
lar disks.  In  fact,  each  x.  .  is  a  b,  ,  of  another  C,  and  one  can 

ij        kh  k 

"walk"  on  the  boundary  of  I  (or  I)  in  this  way.  The  areas  of  I 
and  I  are  oriented  sums  of  the  areas  defined  by  the  arcs  in  V  A 


-10- 


and'XTv.  and  some  set  of  fixed  coordinate  axes  and  thus  can  be 
1   1 

computed  in  O(logn)  time  in  the  CREW  PRAM.  The  data  structure 
constructed  by  the  algorithm  implies  an  0(1)  time  for  answering 
queries  of  the  form  "given  the  point  (x,y)  does  it  belong  to  the 
boundary  of  I  or  the  boundary  of  T."  . 

Note  also  that  the  algorithm  just  presented  is  of  Odoglogn) 
time  in  the  CRCW-PRAM  model,  since  logical  AND  takes  0(1)  time 
there  and  "max"  takes  O(loglogn)  time. 

6.  An  argument  for  optimality  of  the 
CREW-PRAM  algorithms. 

By  using  our  algorithm  we  can  determine  the  area  of  the 
union  of  n  circular  disks  in  O(logn)  time.  Given  a  sequence  of 
n  numbers  in  the  field  of  reals  or  in  any  subfield  of  relative- 
ly long  integers,  call  them  x^,...,x  ,  we  construct  the  disks 
of  radii.       Vx.,...,N/T  and  center  coordinates  (3k+l  maxVx.,0) 
(which,  obviously,  are  non-intersecting).  The  area  of  their 
union  is  then  n- sum ( x ,,..., x  ).  Any  algorithm  to  compute  the 

In         J  a 

area  of  the  union  in  o(iogn)  time  would  also  compute  sumlx  ,...,x 
in  o(iogn)  time.  But  this  contradicts  to  a  result  of  I.  Parberry 
(see  [Pa,  1986]  )  . 


7.  The  number  of  arcs  in  the 
boundary  of  the  intersection. 

Theorem  7 . 1   The  boundary  of  the  intersection  I  of  n  circular 
disks  in  the  plane  consists  of  at  most  2n-2  arcs,  for  any  n^2 


Proof 

We  need  three  helpful  lemmas: 

Lemma  7 . 1   Let  arcs  x,  y   and  y^  be  in  the  boundary  B(I) 
of  the  intersection  I  of  n^2  circles.  Let  x  be  between  y^  and 
y-  and  let  it  intersect  with  y   at  point  A  and  with  y„  at  point 
B.  Let  y   and  y_  belong  to  the  circumference  C(Y)  of  disk  Y  and 
X  to  the  circumference  C(X)  of  disk  X. 
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Then,  C(X)  has  only  one  arc  (i.e.  x)  in  B(I). 
For  proof  of  Lemma  7.1  see  Appendix  A2 

Lemma  7.2   Let  arcs  x  ,  x-  (of  C(X)  of  disk  X)  and  arcs 
y.,  y_  (of  the  circumference  C(Y)  of  disk  Y)  belong  to  the  boun- 
dary B(I)  of  the  intersection  I.  Then,  if  we  start  from  x   and 
order  arcs  on  B(I)  by  the  order  we  meet  them  when  proceeding  on 
B(I)  counterclockwise,  either  both  y  ,  y^  will  appear  before  Xp , 
or  both  will  appear  after  x_. 
For  proof  of  Lemma  7.2  see  Appendix  A3 

Lemma  7.3   For  any  'a>,2  ,    there  exists  a  disk  (out  of  the  n 
disks  participating  in  B(I))  which  has  only  one  arc  on  B(I), 
provided  that  I  does  not  collapse  to  a  point. 
Proof 

Assume,  for  sake  of  contradiction,  that  all  circles  contri- 
bute at  least  two  arcs  each  in  B(I).  Consider  two  arcs,  Xq.,x„p 
of  such  a  circle,  X„ ,  and  let  us  look  at  the  one  of  the  two  con- 
nected pieces  of  B(I)  which  joins  x_   with  x^„  (call  it  A^^). 


01 


'02 


0 


Let  X   be  the  circle  whose  arc,  x..,  connects  to  x^   of  X^ 
and  belongs  on  A  .  Then  (by  Lemma  2)  any  other  arc  x  _  of  X 
must  lie  on  A^  too.  Let  x^p  be  the  arc  of  X   on  A^  such  that  no 
other  arc  of  X   is  between  x    and  x.„  on  A^.  Let  A   be  the  pie- 
ce of  A^  between  x    and  x.„. 


Figure  6 
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I,  ,.  We  distinguish  3  situations: 
k+  1 

Figure  7 


I 

\ 


w 


Case  I 


Case  II 


Case  III 


Case  I   The  number  of  vertices  of  B(I  )  left  out  by  w  is  zero 
In  this  case,  w  increases  the  number  of  arcs  by  2,  since  it 
replaces  a  single  arc  of  a  single  circle  by  3  arcs. 

Case  II   The  number  of  vertices  of  B(I  )  left  out  by  w  is  1. 
In  that  case,  w  increases  the  number  of  arcs  by  1. 

Case  III   The  number  of  vertices  of  B(I  ),  left  out  by  w,  is 

^1.  In  that  case,  w  does  not  increase  the  number  of  arcs  in 

B(I,  ,  )  . 
k+  1 

So,  the  biggest  increase  is  that  of  case  I.  Hence,  the 

number  of  arcs,  P   -,  in  B(I    )  becomes  (in  any  case) 

Pk.1^Pk^2i2k-2.2 

=  2^+1  )-2. 

So,  the  theorem  holds  for  n=k+1,  hence  for  all  n^2. 
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APPENDIX   A1 


Part  a 


Let  e  .  =Right  (b .  . ;E.  )  and 
e.^  =  Right  (b.^^;E.) 

Assume  (for  contradiction)  that  A.  .OA  .  ,  diZ    and  A.,  /A.  ..  This 

ij   Ik         Ik   ij 

means  that 

[b.  .  ,e.  ]  A  [b.  ,  ,e.  -\    ^2 
■■  1 J   ir  -■  '  •■  Ik   is  -■ 

By  the  definition  of  e.   we  must  have 

ir 

b  .  .  :  <  e  .  ,  e  .  *> 

1 J     ir   IS  ' 

Also,  similarly,   b.,  :<e.  ,e.  S 

■'     ik  ^  IS   ir  / 

Since  jep(i)  we  get    b..:<e.  ,b.  ^  Vm;^i,j 

ijimim*        "^ 

Also,  since  kep(i)  we  get   b.,  :^e.  ,b.  >  Vm;^i,k     (EQ8) 

^  ik  ^  im   im ' 

From  (EQ5 )  ,  (EQ6 )  we  get 

(EQ9) 

From  (EQ7)  we  have  b .  . :  ^ e .  ,  , b .  ,  S  (EQIO) 

1 J  ^  Ik  '  ik  ^ 

So,  there  are  two  possibilities  for  the  position  of  e 


(EQ5) 
(EQ6) 
(EQ7) 


b:.:<e.  ,b.,,e.  > 
ij     ir   ik   is'^ 


ik 


Either 


or 


b  .  .  :  ^e  .  ,  ,  e  .  ,  b  .  ,  *> 

1 J  ^  ik   ir   ik  ' 

b..:<e.  ,e.,,b.,S 
ij  ^  ir   Ik   ik' 

The  first  possibility  can  not  happen  because  of  the  defi- 


nition of  e.  . 
ir 


(EQ1  1 


We  conclude  that    b.  .:^e.  ,e.,  ,b   *> 

ij     ir '  ik'  ik  ' 

From  (EQ8)  we  get   b   :<e   ,b  >.By  similar  reasoning  we  con- 

IK       1  J    1  J 

elude  that 

^k=<^s'^ij'^j>  (EQ12) 

But  then,  clearly,  [be]  and  [b.,,e.  ]cannot  intersect  which 

-i-J   ir  IK   IS 

is  a  contradiction 


n 


The  proof  of  Part  2  is  very  similar 


QED 
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Proof  of  Lemma  7 . 1 


APPENDIX   A2 


Figure  A 


C(X) 


B(I) 


Set  I  is  a  convex  set.  Hence,  the  rest  of  the  boundary  of 
I  (except  x)  is  totally  within  circle  X 

□ 


APPENDIX   A3 


Proof  of  Lemma  7.2 


Consider  the  intersection  I(X,Y)  of  disks  X,Y 


Figure  5 


Let  A,B  be  the  points  of  intersection  of  C{X),  C(Y).  The 
line  AB  totally  separates  all  the  points  of  the  arc  AB  of  C(Y) 
from  all  the  points  of  the  arc  AB  of  C(X)  (except,  of  course 
for  A,B).  If  the  arcs  x  ,Xp,y  ,y   were  "intermixed"  in  B(I)  then 
no  straight  line  would  separate  both  x  ,x_  from  both  y  ,y_. 
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